Au4o: A Large Tetrahedral Magic Cluster 
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40 is a magic number for tetrahedral symmetry predicted in both nuclear physics and the electronic 
jellium model. We show that Au40 could be such a a magic cluster from density functional theory- 
based basin hopping for global minimization. The putative global minimum found for Au4o has a 
twisted pyramid structure, reminiscent of the famous tetrahedral Au20, and a sizable HOMO-LUMO 
gap of 0.69 eV, indicating its molecular nature. Analysis of the electronic states reveals that the gap 
is related to shell closings of the metallic electrons in a tetrahedrally distorted effective potential. 

PACS numbers: 61. 46. Be, 36.40.Cs, 31.15.E- 



Gold is a magic clement that constantly brings us sur- 
prises. The relativistic effect has been demonstrated 
to be a dominating factor that distinguishes gold from 
the same-group elements such as copper and silver—. 
Many intriguing structures and properties are found 
for nanometer-sized gold clusters. One example is the 
transition from a two-dimensional structure to a three- 
dimensional one around a size of 10 atoms^. Another 
is the unique catalytic activity of nanometer-sized gold 
clusters^. 

Among all gold nanoclusters, Au2o is probably the 
most famous one^. It has a simple tetrahedral symme- 
try with every gold atom on the surface. The tetrahe- 
dral symmetry is very stable and is the ground state 
both in charged and neutral clusters^. The symmetry 
is still retained if one atom is removed to form Auig^ 
and is proposed to be present in chemically synthesized 
Au2o(PPh3)8 clusters^. 

It is common for gold clusters that the frontier or- 
bitals around Fermi energy are derived from Au(6s) 
atomic orbitals and form delocalized states distributed 
over the whole cluster—. Very similar to the stability 
of noble gas atoms in the periodic system of elements, 
a closed electronic shell built up from these delocalized 
states improves the energetics and gives more chemi- 
cal stability, i.e. reduced reactivit}*^. This effect also 
leads to enhanced stability for specific sizes in protected 
gold clusteraiSiii. Spherical shell closings are especially 
prominent as these are known to lead to large gaps be- 
tween highest occupied (HOMO) and lowest unoccupied 
molecular orbitals (LUMO), the signature of chemical 
stability. 

In the case of deformations of the background po- 
tential, the Jahn- Teller effect can lead to stabilizations 
where spherical shell closings are not available^. A spe- 
cial case here is octupole deformations, where in partic- 
ular tetrahedral deformations produce large gapsi^. 

The appearance of tetrahedral deformations is pro- 
posed in nuclear physic&i^ where the delocalized particle 
picture in an effective background potential had origi- 
nated before it was adopted in cluster physics. While the 
experimental observation of the tetrahedral deformation 
in nuclei seems to be under debat o^^'^^ , the importance 



of tetrahedral symmetry in cluster physics in the case of 
Au2o is beyond questior>i^. 

Magic tetrahedral metal clusters were proposed in 
the jellium modeU^. In particular, in the essen- 
tially parameter-free ultimate jellium model, where the 
smeared out nuclear density exactly follows the electronic 
density, the 40 electron ground state shows a tetrahedral 
deformatiortii. In spherical symmetry the electrons' an- 
gular momentum is conserved and the relative energies 
of different angular momentum shells (and also the gaps 
between them) depend on the effective radial potential. 
There is a large gap for a harmonic radial potential at 40 
electrons. This gap decreases when the potential becomes 
more box like, however—. Here a tetrahedral deformation 
can open the gap again. To our knowledge, tetrahedral 
metal clusters larger than the rather trivial case of four 
atoms have been observed only in the example of Au2o. 
In this Letter, we show from density functional theory- 
based global minimization that the larger Au4o has a 
twisted trigonal pyramid structure of quasi-tctrahedral 
symmetry, the first case of a metal cluster of tetrahedral 
shape beyond Au2o. 

Although gold clusters with 20 atoms or less have been 
extensively studied, we know relatively little of the struc- 
tures of larger gold clusters. Au clusters with 30 to 60 
atoms would be the key to understanding the transition 
from the molecular behavior of a small cluster to the 
metallic bulk. The clusters Au^2^'^ and AU34— have 
been proposed to have a core-shell structure, instead of 
being hollow or planar. More recently, the global minima 
of Au28 to Au35 were explored in comparison with exper- 
imental photoelectron spectrap^i. The authors found that 
the global minima are amorphous in nature with an Au4 
tetrahedron core and a much bigger outer-shell for Aun 
with TV > 32. Using an empirical potential for global 
minimum search followed by density functional theory 
(DFT) calculations, Garzon ct al. found that the most 
stable structure of Auag is of C's symmetry with an Au5 
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This structure was found to be slightly lower 
in energy than the high-symmetry truncated octahedron 
(by about 0.3 to 0.6 eV, depending on the choice of DFT 
functional^i^^). Both the Cg and octahedral structures 
are metallic (that is, their HOMOs are not completely 
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filled). Using a strategy similar to Garzon et al.'s, Tran 
and Johnston found a structure for AU40 with a distorted 
truncated octahedroii^i. 




Basin-hopping steps 

FIG. 1. Energy landscape of a DFT-based basin-hopping 
search for the global minimum of Ausg, started with the oc- 
tahedral structure. 

Puzzled by the metallic nature of Auas of the state-of- 
the-art models and expecting a high-symmetry structure 
for Au4o, we set out to find the global minima of Auag 
and Au4o, hypothesizing that Au4o's structure will build 
upon that of Auag. Instead of using the empirical poten- 
tials to search for global minima before DFT optimiza- 
tion, as previously done for Auag and AU40 , we use DFT 
geometry optimization at the GGA-PBE level^ directly 
in our basin-hopping global-minimum searcb^^. This ap- 
proach has been quite powerful for exploring the energy 
landscape of nanoclusters^Ii^. What distinguishes our 
work from previous DFT-based basin-hopping search for 
the similar-sized gold clusters is that we run the basin- 
hopping procedure for much more steps (over 1000) to 
explore the energy landscape. 

Fig. 1 shows our DFT-based basin-hopping search for 
the global minimum of Auag, started with the highly sym- 
metric octahedral structure. One can see that the octahe- 
dral structure (Fig. la) was transformed into much less 
symmetric configurations of lower energy. The energy 
lowering is quite substantial, more than 2 eV, instead of 
the 0.3 to 0.6 eV lowering found by Garzon et al. for 
their Cs mode l^^'^'^ . After a local minimum with one 
adatom sticking out (Fig. lb), we found two putative 
global minima (Fig. Ic and Id) which are almost degen- 
erate in energy (within 30 meV). The two structures have 
a similar construction: an Au4 core and an AU32 shell, 
with two Au adatoms sticking out; but one with Ci sym- 
metry (Fig. Ic), the other with C2 symmetry (Fig. Id). 
What is unique about the two structures is that they 
both have a sizable HOMO-LUMO gap. indicating their 



moleculc-like stability in the gas phase. The Ci structure 
has a gap of 0.66 eV and the C2 structure 0.84 eV, un- 
like the metallic nature of the octahedral structure and 
Garzon et al.'s Cs mode P^-^'^ . Hence we confirmed that 
the nanometer-sized Auag still behaves as a molecule. 
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Basin-hopping steps 

FIG. 2. Energy landscape of a DFT-based basin-hopping 
search for the global minimum of Au4o. 




FIG. 3. The global minimum for Au4o featuring a twisted 
pyramid with a missing corner: (a) top view, the missing 
corner is at the top of the figure; (b) side view, the Au4 tetra- 
hedral core is highlighted in red. 

The two Au38 minima are also interesting in that they 
have two adatoms on the Au4@Au32 core shell. This 
construction shares some similarity with the AU35 struc- 
ture where one adatom is on the Au4@Au3o core shell^i. 
What Au38's structure suggested to us is that Au4o's 
structure can build upon the Au4@Au32 core-shell frame- 
work but with four adatoms placed in tetrahedral sym- 
metry. This idea led us to propose an initial guess for 
Au4o (Fig. 2a) with two more adatoms manually added 
to one of the Ausg's two candidate structures (Fig. Ic). 
Started with this initial guess, we performed DFT-based 
basin-hopping search for over 1000 steps. Interestingly, 
twisted pyramid structures evolved out. The first one 
evolved out is a twisted trigonal pyramid with a missing 
corner (Fig. 2b and Fig. 3); the core is still an Au4 tetra- 
hedron (Fig. 3b); overall, the cluster has Ci symmetry. 
The second one evolved out is a twisted trigonal pyramid 
with a missing core atom (Fig. 2c and Fig. 4); this struc- 
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ture has C3 symmetry, with an Aus triangle core (Fig. 
4b). The C3 symmetry can also be clearly seen from the 
base of the pyramid (Fig. 4c). The top three layers of 
the pyramid share the same substructure as the famous 
tetrahedral Au2o. Both Au4o structures are chiral as was 
first proposed for Aus^^. 




FIG. 4. An isomer of Au4o featuring a twisted pyramid with 
C3 symmetry: (a) top view; (b) side view, showing the Au3 
triangle core; (c) bottom view, showing the base of the pyra- 
mid. 

To gain a deeper insight into the consequences of 
the tetrahedral symmetry we now discuss the electronic 
structure and energetics of the two Au4o isomers. We 
found that the Ci isomer is the most stable among all 
the configurations explored; it has a HOMO-LUMO gap 
of 0.69 eV. The C3 isomer is only 0.15 eV higher in energy 
and has an even larger gap of 0.85 eV. So both isomers 
are stable molecules in nature, given these rather large 
gaps in clusters of this size. Moreover, we found that 
if one manually moves the apex atom of the Ci isomer 
to the missing corner (Fig. 3a), the resultant structure 
is only slightly higher in energy (by 0.07 cV), indicating 
the robustness of Au4o's tetrahedral shape despite the 
multiple isomers close in energy. 

The large HOMO-LUMO gaps of AU40 are related to 
tetrahedral symmetry. In perfectly spherical clusters one 
finds each DFT Kohn-Sham orbital to be in an unique 
angular momentum eigenstate relative to the cluster's 
center of mass. Due to the deformation of the nuclear 
background away from spherical symmetry these states 
are not clean anymore. A tetrahedral deformation as 
present here can be described by an effective potential of 
the forroi^i^a 

V{V) = V{r) [1 + a32(T3,+2 - r3,_2))] (1) 

where r = |r|, the 13.^2 a-i'c spherical tensor operators 
and 0132 is a constant describing the degree of the defor- 
mation. The exact form of the T3,±2 is not important 
for our purpose; one only has to note that these opera- 
tors couple angular momentum eigenstates with angular 
momentum projections that differ by ±2 exclusively. 

With these considerations wc analyze the electronic 
structure of the tetrahedral gold clusters in Fig. O^. 
We project the DFT Kohn-Sham orbitals onto spheri- 
cal angular momentum eigenstates relative to the clus- 
ter center of mass to extract their delocalized nature^^. 
The usual nomenclature characterizing the states' angu- 
lar momentum is similar to atomic physics S, P, D, . . . , 
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FIG. 5. (Color online) The density of states relative to Fermi 
energy projected on angular momentum eigenfunctions rela- 
tive to the clusters center of mass. In a, b) only the Au(6s) 
electrons are treated as valence electrons, whereas in c, d) also 
the Au(5d) electrons are allowed to rearrange. 



where capital letters distinguish from atom centered an- 
gular momenta. The principle quantum number n gives 
the energetic ordering which is connected to the (n — 1) 
number of radial nodes of the corresponding eigenfunc- 
tions. In spherical clusters one would expect the delocal- 
ized states to fill lS2lp6lDi"2S2 orbitals for 20 electrons 
and additionally the 1F^^2P^ orbitals for 40 electrons. 

To obtain a clean picture of delocalized states we first 
fix all the electrons except the Au(6s) electrons in a frozen 
core approximation. The projected density of states 
(PDOS) for Au2o obtained in this way is presented in 
Fig. Oi. The states of lowest energy are the delocalized 
IS and three IP states, similar as in the spherical case. 
This is understandable as the magnetic quantum number 
of these states is 0, ±1 and hence these could only cou- 
ple to higher angular momentum states far away in en- 
ergy through the T3_±2 operators in eq. H]). Then there 
follows a block of 4 states with mixed S/P/D/F sym- 
metry due to the T3^±2 coupling, immediately followed 
by two states of pure D symmetry that form the Au2o 
HOMO. After the substantial gap of 2.06 eV, the Au2o 
LUMO consists of 3 states with S/P/G symmetry and is 
followed by a block of 3+2 states of dominant F symme- 
try. After these there is another large gap. This analysis 
clearly shows how sparse the delocalized electronic states 
are distributed under tetrahedral deformations and that 
the electronic system can make gain from closing shells 
in the corresponding symmetry. 

Analyzing the Au4o states in the same way leads to a 
very comparable picture as shown in Fig. [SJj. The rela- 
tive gaps between the blocks of states gets smaller due 
to the larger size of the cluster, but the symmetry of the 
states is similar. In this cluster the rather large HOMO- 
LUMO gap is between the F/P symmetry dominated oc- 
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cupicd states and the G symmetry dominated unoccupied 
states. Including Au(5d) states as valence electrons into 
the calculation does not change the picture of the frontier 
orbitals around the Fermi energy as shown in Fig.[5j; and 
d. A comparison with the C3 isomer that has a clearer 
tetrahedral structure indicates that in terms of symmetry 
the HOMO of the Ci isomer belongs rather to the LUMO 
block of states with dominant P/G symmetry. As a con- 
sequence the gap of the C3 isomer is even larger than that 
of the Ci isomer. Finally we have analyzed the deforma- 
tion of the s-valence electron density as it was done for 
the near tetrahedral shape of Na4o^. While we obtain 
for the largest distortion parameter S3 = 0.04 for Na4cr^ 
in agreement with ref. |3l|, ^3 = 0.26, 0.33 for the Ci, C3 
isomers of Au4o respectively, accounting for the much 
larger tetrahedral distortion present in the gold clusters. 

In summary, we found AU40 to be a magic cluster with 
a quasi-tetrahedral symmetry. It has a twisted pyra- 
mid structure discovered from DFT-based basin hopping 
for global minimum search and built upon the putative 



global minima of Ausg . This cluster is a manifestation of 
the enhanced stability due to the tetrahedral symmetry, 
predicted both in nuclear structure and by the jellium 
model. Analysis of the delocalized electrons in Au4o con- 
firms the shell-closing picture by the tetrahedral symme- 
try, similar to that of Au2o. The delocalized 6s electrons 
and the complex energy landscape for clusters such as 
Au4o cannot be accurately described by empirical poten- 
tials, thereby making DFT-based global-minimum search 
a necessity. 
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